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[i] A model combining low-frequency complex conductivity and high-frequency 
permittivity is developed in the frequency range from 1 mHz to 1 GHz. The low-frequency 
conductivity depends on pore water and surface conductivities. Surface conductivity is 
controlled by the electrical diffuse layer, the outer component of the electrical double layer 
coating the surface of the minerals. The frequency dependence of the effective quadrature 
conductivity shows three domains. Below a critical frequency f p , which depends on the 
dynamic pore throat size A, the quadrature conductivity is frequency dependent. Between f p 
and a second critical frequency f d , the quadrature conductivity is generally well described 
by a plateau when clay minerals are present in the material. Clay-free porous materials with 
a narrow grain size distribution are described by a Cole-Cole model. The characteristic 
frequency f d controls the transition between double layer polarization and the effect of the 
high-frequency permittivity of the material. The Maxwell-Wagner polarization is found to 
be relatively negligible. For a broad range of frequencies below 1 MHz, the effective 
permittivity exhibits a strong dependence with the cation exchange capacity and the specific 
surface area. At high frequency, above the critical frequency f d , the effective permittivity 
reaches a high-frequency asymptotic limit that is controlled by the two Archie's exponents 
m and n like the low-frequency electrical conductivity. The unified model is compared with 
various data sets from the literature and is able to explain fairly well a broad number of 
observations with a very small number of textural and electrochemical parameters. It could 
be therefore used to interpret induced polarization, induction-based electromagnetic 
methods, and ground penetrating radar data to characterize the vadose zone. 

Citation: Revil, A. (2013), Effective conductivity and permittivity of unsaturated porous materials in the frequency range 1 mHz- 
1GHz, Water Resour. Res., 49, doi: 10. 1029/20 12WR0 12700. 



1. Introduction 

[2] There are a number of electrical and electromagnetic 
(EM) geophysical methods providing nonintrusively infor- 
mation relative to the texture, water content, and interfacial 
electrochemistry in the shallow subsurface in an extremely 
broad range of frequencies (1 mHz-1 GHz). The relevant 
arsenal of nonintrusive techniques include self-potential, 
direct current (DC) resistivity, frequency and time-domain 
induced polarization, frequency and time-domain induction 
EM methods, ground penetrating radar (GPR), and the seis- 
moelectric/electroseismic methods. These nonintrusive 
methods have many applications in hydrogeophysics to 
assess ground water resources and to monitor ground water 
quality or for vadose zone hydrogeology and agriculture 
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just to cite few examples [e.g., Linde et al., 2006]. That 
said, there is not a single model able to explain the complex 
conductivity of unsaturated materials as well as their effec- 
tive permittivity and their frequency dependences in the 
broad range of frequencies of interest in geophysics 
(1 mHz-1 GHz). 

[3] In the low-frequency domain (<10 kHz), new mod- 
els have been developed recently to describe the complex 
conductivity of saturated and unsaturated sands and clayey 
materials [Revil and Florsch, 2010; Revil, 2012; Revil 
et al., 2012a, 2012b]. There is a vast literature regarding 
this topic in interfacial electrochemistry and colloidal 
chemistry [e.g., Schwarz 1962; Dukhin and Shilov, 2002; 
Grosse, 2002] as well as in geophysics [Tarasov and Titov, 
2007]. At very high frequencies, there are also predictive 
models describing the high-frequency permittivity of po- 
rous rocks. These models are empirically based like the 
Topp model for soils [Topp and Reynolds, 1998], based on 
upscaling methods like the volume-averaging method 
[Pride, 1994], or based on the effective differential medium 
approach or mixing laws (see discussion in Jones and Or 
[2002], Cosenza et al. [2003], Robinson and Friedman 
[2003], and Miyamoto et al. [2005]). Note that with the vol- 
ume-averaging approach, the different phases are consid- 
ered symmetrically while with the differential effective 
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medium approach, one phase has to be chosen as the host 
(for instance, grains immersed in a background fluid or 
fluid inclusion in a matrix). 

[4] The situation is much more confused in the interme- 
diate range of frequencies (10 kHz-10 MHz). Typically, a 
number of authors in geophysics have extended the high- 
frequency behavior to lower frequencies by adding a con- 
tribution related to the Maxwell- Wagner polarization 
mechanism [Chen and Or, 2006]. Other approaches have 
relied on using the so-called standard model of colloidal 
chemistry (developed initially by Dukhin and Shilov and 
based on the polarization of the diffuse layer, see Grosse 
[2002]) and applying this model to low-frequency data sets 
[Garrouch and Sharma, 1994]. To fit the data, a stretch of 
the textural parameters is usually required and the obtained 
set of interfacial properties are not compatible with the tri- 
ple layer complexation models commonly used in surface 
electrochemistry of minerals. 

[5] In this paper, following my previous works [Revil 
and Florsch, 2010; Revil, 2012; Revil et al., 2012a, 
2012b], I claim that the diffuse layer is not the main con- 
tributor to low-frequency polarization of clayey materials. I 
also claim that the Maxwell- Wagner polarization is not ei- 
ther the dominating mechanism of polarization at interme- 
diate frequencies. An extremely simple model is developed 
based on the low-frequency polarization of the Stern layer 
(the inner layer of the electrical double layer). This model 
is proposed for saturated and unsaturated clayey materials 
(water being the wetting phase). It explains for the first 
time an important amount of literature data in the fre- 
quency range 1 mHz to 1 GHz in a very simple way. 

2. Theory 

[6] I consider an unsaturated porous material as shown 
in Figure la with an electrical double layer coating the sur- 
face of the grains (Figure lb). In the classical literature in 
geophysics, it is customary to see a polarization sketch like 
the one shown in Figure 2a [e.g., Hasted, 1961]. In the fre- 
quency range 1 mHz-1 GHz, the apparent dielectric con- 
stant of a porous material is explained by two dissipation 
mechanisms only: the Maxwell Wagner polarization 
(called the j3 polarization in electrochemistry, see Grosse 
[2002]) and the polarization of the water molecules past the 
gigahertz frequency (called the 7 polarization). That said, 
we know that in induced polarization, the polarization of 
the electrical double layer plays an important role below 
10 kHz. This nondielectric dispersion mechanism is associ- 
ated with the electromigration and accumulation/depletion 
of charge carriers at discontinuities in the migration path- 
ways of the charge carriers. This mechanism is called the a 
polarization in electrochemistry, and its role is such that it 
may be the dominating mechanism of polarization for a 
very broad frequency range as discussed in this paper (see 
Figure 2b). 

2.1. Effective Conductivity and Permittivity : 
Definitions 

[7] I first define the concepts of effective conductivity, 
effective permittivity, and effective quadrature conductiv- 
ity. Ampere law is given by [e.g., Vinegar and Waxman, 
1984], 
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Figure 1. The REV. (a) Sketch of the REV comprised 
between two reservoirs, (b) Sketch of the electrical double 
layer at the solid water interface. The double layer includes 
the Stern layer and the diffuse layer. The charges of these 
two layers compensate the charge on the mineral surface 
(o plane). The d plane corresponds to the separation plane 
between the Stern and the diffuse layers, (c) Sketch of the 
solid-water and water-air interfaces. 



V x H = J + 



9D 

~dt' 



(1) 



where J = a* E denotes the conduction current density (in 
A m -2 ), a* = o 7 + ia" represents the complex conductiv- 
ity, 3d = dD/dt corresponds to the displacement current 
density (A m~ 2 ), D = eE represents the dielectric displace- 
ment (in C m~ 2 ), H denotes the auxiliary magnetic field 
(A m _1 ), and e denotes the permittivity or dielectric con- 
stant (in F m _1 ) of the material. I consider an harmonic 
external electrical field E = E 0 exp(— iojt) (i = v^— T is the 
pure imaginary number, / is the frequency in Hz, w = 2nf 
the angular frequency (in rad s _1 ), and E 0 denotes a con- 
stant amplitude of the alternating electrical field). The total 
current density J, is the sum of J and Jj. This yields J t = 
(a* - iujs)E [Vinegar and Waxman, 1984]. This total cur- 
rent density can be written as, 



J, 



(2) 



where er* ff = a e ff — ;a>£ e ff is the effective complex conduc- 
tivity and <7 ef f and e e ff are real positive frequency-depend- 
ent scalars defined by 



<7eff = o , 
£eff = e' - a"/u. 



(3) 
(4) 



It is also possible to define an effective quadrature conduc- 
tivity (in S m -1 ) that is directly related to the effective per- 
mittivity as, 



°eff = W£ eff 



-ct" + ojs'. 



(5) 
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a. Traditional Sketch of dielectric dispersion 
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Figure 2. Dielectric dispersion for polarization in 
charged porous media, (a) Classical representation from 
Poley et al. [1978] and found in a number of papers in geo- 
physics, (b) Representation of the three types of polariza- 
tion discussed in the text. The a polarization corresponds 
to the polarization of the electrical double layer, the (3 
polarization corresponds to the Maxwell-Wagner polariza- 
tion, and the 7 polarization concerns the polarization of the 
water molecules. 

Note that in the conventions used above, <r" ff is positively 
defined while a" is negatively defined. Two phase angles 
can be defined. The former is determined from the complex 
conductivity tan</> = a" /a 1 and the second from the effec- 
tive parameters tan6 = a" ff /a e ff = u;e e ff /<7 e ff (note that 
<f> < 0 while 6>0). 

2.2. Conductivity and Permittivity : 
A Volume-Averaging Approach 

[8] Pride [1994] developed a volume-averaging 
approach used to upscale the Nernst-Planck equation to the 
scale of a representative elementary volume (REV) of a 
porous rock. He was however concerned only by the satu- 
rated case and motivated by the modeling of the seismoelec- 
tric effect. We use this approach to upscale the Nernst- 
Planck equation to a porous material saturated by two immis- 
cible fluid phases (water and air). The local current density 
j w (in A m -2 ) occurring in the water phase is determined 
from the Nernst-Planck equation, which is written, in ab- 
sence of macroscopic concentration gradients, as 

N 



where e w denotes the local electrical field (in V m _1 ), yV 
denote the number of ionic species, q { the charge of species 
i, and C, their concentration. The mobilities 6, (in N s m _1 ) 
entering equation (6) are related to the classical ionic mobi- 
lities fit (in m 2 s -1 V -1 ) by b t — and to the ionic 
self-diffusion coefficients Z), by Z), = k^Tbi, where kt 
denotes the Boltzmann constant (1.3807 x 10" 23 J K" 1 ), 
and T the absolute temperature (in K). There is also a cur- 
rent density associated with the grains, j s due to the electri- 
cal double layer coating the surface of the insulating grains. 
The boundary conditions at the interface S sw separating the 
solid phase and the pore space are 

n. v ■ (e, - e w ) = 0, (7) 

■va-w = o, (8) 

where n s denotes the normal unit vector to S sw directed 
from the pore to the solid phase (Figure lc). The pore space 
is filled with two immiscible fluids, the water and air 
phases (water is assumed to be the wetting phase). The 
boundary conditions at the interface S aw separating the non- 
wetting phase (air) and water are 

n a ■ (e a - e w ) = 0, (9) 
n«-(i a -y=0, (10) 

where n fl denotes the normal unit vector to S aw directed 
from the water to the nonwetting phase (Figure lc). In the 
following, I will consider that the surface charge density at 
the interface between water and the nonwetting phase can 
be neglected. The reason to neglect this contribution is 
because there is still a gap of knowledge in describing the 
conductivity and polarization of this interface, which is, 
however, described by an electrical double layer [Leroy 
et al., 2012] and therefore likely characterized by surface 
conduction and polarization. 

[9] The microscopic equations given above are now 
averaged at the scale of an REV. The REV corresponds to 
the rock volume located in between two large-parallel 
circular disks of area A separated by the distance L. This 
is typically the case of a jacketed cylindrical sample in 
the laboratory. We assume that there is a macroscopic 
electrical potential difference Aip. When this potential 
difference is divided by L, one obtains the equivalent 
macroscopic field perpendicular to the end faces of the 
REV. We note z the unit vector normal to the end faces 
(Figure la). For example, for the electrical field at the 
pore or grain scale, we have V x e^ = 0, so the electrical 
fields e^ can be derived from scalar electrical potentials 
e^ = — Vf/^ where the index £ defined the phase (w for 
water, s for solid, a for air). The macroscopic electrical 
field is written as 

z-E = --^, (11) 

where Aip = ip^L) - ^(0) denotes the difference of elec- 
trical potential between the two reservoirs (Figure la). In 
the pore water phase, the fundamental Laplace problem is 
defined by 
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V 2 r w = 0, in v w , 

n s ■ Vr w = 0, on 5" sw , 

n fl • Vr„, = 0, on S^, 

L, on z = L 
0, on z — 0, 



(12) 
(13) 
(14) 

(15) 



where S sw and S aw denote the solid-water and air-water 
interfaces, respectively (Figure lc), and V w denotes the 
water phase volume. The electrical potential in the water 
phase can be written as ip w = T w Aip/L and the local elec- 
trical field by e vv = —VT w Aip/L. A similar boundary-value 
problem of the normalized effective potential T a can be 
defined in the nonwetting phase (V a denotes the air phase), 



V 2 r„ = 0, in V a , 

n a • Vr a = 0, in S aw , 

L, on z = L 
0, on z = 0. 



(16) 
(17) 

(18) 



[10] The volume average of a quantity a f (a^ is a scalar, 
a vector, or possibly a tensor) is 



(19) 



where is the volume of the £ phase within the REV. Slat- 
tery's classical theorem states [Slattery, 1967], 



<Va ? ) = V<a e )+^ J n^dS, 



(20) 



(no summation on the indices) where n ?(T represents the 
unit outwardly directed normal vector for the £ phase and 
Sfr represents the interfacial area contained within the 
averaging volume. Applying this to the water, solid, and air 
phases and using the unit vectors defined in Figure lc, I 
obtain, 

(Va*) = V(a„,) + yj n s a w dS + yj n a a w dS, (21) 

s, s w 

(Va,) = V(a s ) - iy n s z s dS, (22) 

s, 

(Va a ) = V(a a ) - iy n a a a dS, (23) 



where dS is an infinitesimal surface volume element. The 
volumetric phase average a^ and the total volumetric aver- 
age A are rerelated to each other by, 



a? = (»£)//?> 



(24) 
(25) 



where fc denotes the volumetric fraction of phase £(fjr = <j>, 
(f> denotes the connected porosity, and / s = 1 — (f> is the so- 
lidity, f w = s w (f> = 9 is the volume fraction of water, also 
called the water content in hydrogeology, and f a = 
(1 — s w )<j> denotes the volume fraction of the nonwetting 
phase). Therefore, we have, 



A = (1 - </>)n s + s w 4>a w + (1 -s w )(j>A a 



(26) 



[n] We first define the electrical conductivity of the 
pore water, which is given by, 



a w = J2q^b i C i = J2\q i \P l C i . 



(27) 



[12] From equation (26), the average conduction current 
density is given by 



J = -a s J e s dV + e w dV, 



(28) 



v, v w 
J = (1 - <(>)a s e s + s w <f>a w e w . (29) 

[13] The tortuosity of the water phase a w is defined by 

1 



(30) 



— = 1+^- [z-n s T w dS + ±- fz n a T w dS, (31) 



where E = —zAtp/H denotes the macroscopic electrical 
field. From equations (25) and (26), the macroscopic elec- 
trical field is also given by 

E = (e s ) + (e w ) + (e a ), (32) 
E = (1 - cf>)e s + <t>s w e w + 4>(\ - s w )e a . (33) 

The phase average of the electrical field e a can be also 
related to the macroscopic electrical field via the tortuosity 
of the nonwetting phase a a , 



1 - 
e a = — E, 

O-a 



Ota 1 V, 



y j z • n a T a dS. 



(34) 
(35) 



I need now to find a relationship between the tortuosity of 
the pore space a (related to the definition of the electrical 
formation factor F at saturation by F = a /(f)) and the tor- 
tuosities of the wetting and nonwetting phases. If we note 
ey the local electrical field of the fluid phase (comprising 
the wetting plus nonwetting fluid phases), the phase aver- 
age of tf is related to the phase average of the electrical 
field in the wetting phase e w and to the phase average of 
the electrical field in the nonwetting phase e a by 



Cy s w G w + (1 s w )e a 



(36) 
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The electric field tf is related to the macroscopic field E by, 



e f = -E. 

a 



(37) 



Combining equations (30), (34), (A36), and (37), I obtain 
after few algebraic manipulations the following relation- 
ship between the tortuosity of the pore space, the tortuosiy 
of the water phase, and the tortuosity of the air phase, 



1 S W , 1 

a a w a a 



(38) 



Therefore, the tortuosity of the pore space is equal to the 
harmonic average of the tortuosities of each phase weighted 
by their relative saturation. If we multiply equation (38) by 
the porosity, we obtain a relationship between the tortuosity 
of the air phase, the formation factor, and the tortuosity of 
the water phase, 



<ft(l ~ s w) _ _j_ _ <fcw 

a„ F a w 



(39) 



The electrical formation factor can in turn be related to the 
porosity by a power-law function [Archie, 1942], 



(40) 



called the first Archie's law, where m is called the cementa- 
tion exponent (m > 1 is a strict bound and (1 < m < 3 for 
sandstones and sands according to experimental data, see 
Revil et al. [1998, Figure 5]). The conduction current den- 
sity can be written as 



:ff,[E- 



..).S' u C n 



(1 -s w )e a ] +s w <j)a w e w , 



a s E + (j>s w (a w - a s )e w - (1 - s w )tf>a s e a 



- . ft>, 

&s H (o» 



ffj tr. 



(41) 
(42) 

(43) 



The effective electrical conductivity of the material (a) (in 
S m -1 ) is defined by the macroscopic Ohm's law, 



J c = ME. 

Combining equations (43) and (44) yields 



(a) = a s H (a w - a s ) a s . 

Otrti Otn 



(44) 



(45) 



Note that in the case where there is no surface conductiv- 
ity, we can connect equation (45) to the second Archie's 
law by 



Jim (a) = a w = 

<r,^0 a w r 



(46) 



where n is called the second Archie's exponent in the litera- 
ture. In the following, this second Archie's exponent in the 



expression of the electrical conductivity is introduced by 
using the following change of variables, 



ft 



1," 
F s w 



(47) 



From equation (39), I can also use the following change of 
variables, 



(1 



1 



Oia 



(i-O- 



(48) 



Using equations (47) and (48) into equation (45), I obtained 
the following expression for the electrical conductivity of 
the porous medium, 

{o) = *. + jZ&w - <r.) ~ j (1 - (49) 
(a) = -sla w + a s - -a s , (50) 
{<r)=^[s n w a w + (F-l)a s ]. (51) 



In the following, I remove the bar on the pore water and 
surface conductivity terms and drop the average symbol for 
the electrical conductivity to keep the notations as simple 
as possible. The same chain of algebraic manipulations 
could be applied also to complex conductivity and there- 
fore I can write the following expression of the complex 
conductivity of the clayey material, a*, as a function of the 
conductivity of the pore water a w and the complex conduc- 
tivity of the solid phase <r* s = <r^ + ia" s (complex because 
of the polarization the electrical double-layer coating the 
grains), 



(52) 



[14] Note that the pore fluid conduction and surface con- 
duction add in parallel in this model, and therefore the elec- 
trical conductivity of the porous material is a linear 
function of the conductivity of the pore fluid. This is 
expected as the volume-averaging approach treats the solid 
and fluid phases symmetrically, and they are both assumed 
to be continuous through the REV. As discussed in Appen- 
dix A, the situation is different when a differential effective 
medium approach is used. When we are dealing with gran- 
ular media, the solid phase is discontinuous and grains are 
immersed in the background pore water. In this case, the 
differential effective medium theory yields a nonlinear 
relationship between the electrical conductivity of the po- 
rous material and the conductivity of the pore fluid (see for 
instance, Revil et al. [1998] and Revil [1999]). Sands and 
sandstones are somewhere in between these two types of 
models as the diffuse layer is continuous (overlapping from 
grain-to-grain) while the Stern layer is discontinuous. In 
addition, it follows that the diffuse layer does not polarize 
while the Stern layer polarizes and that the Stern layer can- 
not contribute to the DC electrical conductivity. 

[15] Using the displacement current density rather than 
the electromigration current density, an equation similar to 
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equation (52) can be obtained for the high-frequency dielec- 
tric constant using the same chain of algebraic manipulations, 



1 



[& w +(l-£)e tt + {F-l)e s ]. 



(53) 



[i6] Now coming back to equation (52), I need to 
express the low-frequency conductivity of the solid phase 
corresponding to an insulating solid phase coated by a con- 
ductive and polarizable electrical double layer as a function 
of the surface charge density, the mobility of the counter- 
ions, and the tortuosity of the electromigration pathways 
along the surface of the grains. I assume that at low fre- 
quency, only the conduction of the electrical diffuse layer 
participates to the surface conductivity. Indeed, the grains 
are in contact to each other (a very important distinction 
with colloidal suspensions) and the diffuse layer is always 
above a percolation threshold. This assumption has been 
checked experimentally by Vaudelet et al. [201 la, 201 lb]. 

[n] The electrical conductivity of the grains is given by 
the mobility of the counterions (corrected by the tortuosity 
around the grains) times the concentration of the charge 
carriers times their charge. This yields, 



o _ f 3 (+) n S 
a» Vs 



(54) 



where S denotes the surface area of the grains (in m 2 ), Q d is 
the surface charge density of the diffuse layer (in C m -2 ), 
and V s is the volume of the grains (in m 3 ). Note that the mo- 
bility of the charge carriers is corrected for the tortuosity of 
the electromigration of the counterions along the surface of 
the grains in the water phase (the counterions are not allowed 
to pass through the grains that are insulating). The fraction of 
the volumetric charge density associated with the diffuse 
layer is written as 



(1 -fhi)Qv = p s 



(1 -fu)CEC, 



(55) 



where f M (dimensionless) denotes the fraction of counter- 
ions in the Stern layer (the counterions compensating the 
negative charge of the mineral surface are distributed 
between the diffuse layer and the Stern later, see Revil 
[2012] and references therein). Consequently, (1 — fu) 
denotes the fraction of counterions in the diffuse layer. The 
volumetric charge density Q v (in C m -3 ) denotes the total 
charge density of the double-layer per unit volume of the 
pore water (including the Stern and diffuse layers taken to- 
gether). The relationship between Qy and its value at satu- 
ration (Qf, at s w = 1) is Qv = Qy/sw We also have 



0--fu)Qr = Qdw 



where equation (55) has been used. Therefore, the real (in 
phase) low-frequency component of the solid phase con- 
ductivity is given by, 



4 = (l/F4>)s"; /3 {+)Ps (\ -f M )CEC, 



(59) 



where I used the following change of variables 
(l/a w ) <==> (l/F^)s" w resulting from equation (47). I 
replace below the exponent (n — 1) in the expression of the 
complex surface conductivity by the exponent p because 
the true tortuosity that should be used is the one associated 
with pathways along the grains, not strictly speaking, the 
tortosity of the water phase adjacent to the grains (I will 
show however in section 7 that the approximation p m n — 
1 is a good one). This approach is consistent with the one 
derived by Vinegar and Waxman [1984] and Slater and 
Glaser [2003]. The low-frequency (in phase) conductivity 
is given by 



(T 0 : 



F- 1 
FS 



/3 (+) P s (l -/m)CEC 



(60) 



[is] A similar approach can be followed for the high-fre- 
quency electrical conductivity, for which now both the 
Stern and diffuse layers contribute to the overall surface 
conductivity [Revil, 2012; Revil et al., 2012a]. Following 
the same approach as used to derive equation (60), this 
yields 



F(j> 



Ps [/?(+) (l-/j/)+/3f +) /i/]cEC }. 



(61) 



where f3 s ^ denotes the mobility of the counterions in the 
Stern layer (as expected /3? + ) <C /?/ + ), see Revil [2012] for 
clays and Revil et al. [2012a] argued for (3? + ^ = /3( + ) in silica). 

[19] An important parameter in the description of fre- 
quency-domain induced polarization, but more importantly 
in the description of time domain induced polarization, is 
the normalized chargeability. In a time-domain induced 
polarization experiment, the chargeability M characterizes 
how the decay of the secondary voltages when the primary 
current is shut down. This chargeability can be normalized 
to obtain the normalized chargeability M„ = Ma x = cr^ 
-(To, and therefore from equations (60) and (61), the nor- 
malized chargeability is given by 



M n = sP w (J— J^-J Ps^/mCEC . (62) 

We can express also the conductivity terms in term of the 
volumetric charge density Qy, which yields, 



F — 1 



(56) 



Therefore, the surface charge density can be expressed in 
terms of the CEC by 



Q d S={l-f M )Q v V w , 

&|- = p 5 (l-/j,)CEC, 
Vs 



(57) 
(58) 



a 0 = -s w a w + a s 



Too = -«/» + a s 



(63) 
(64) 



in which we have entered the term (F — 1 ) /F ss 1 in 
the expressions of the surface conductivities. In these 
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83 
S 



a. Type A - Berea sandstone. Data from Lesmes and Frye [2001] 



(ay = -o n +(aP) 



Jf 2m 
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b. Type B - Shaly sandstone - Hancock et al, [2012] (Sample #436) 
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Figure 3. Effective quadrature conductivity versus fre- 
quency, (a) Type A. Berea sandstone (data from Lesmes 
and Frye [2001], pH 8, 10" 1 M NaCl solution). The data 
are characterized by three components : a minimum critical 
frequency from which the magnitude of the quadrature con- 
ductivity is going down when the frequency decreases, a 
plateau (over at least four decades in frequency), and a 
high-frequency domain in which there is an increase of the 
magnitude of the apparent quadrature conductivity. Filled 
circles : measurements with four electrodes, filled squares, 
measurements with two electrodes, (b) Type B. shaly con- 
glomerate (this study, pH 7, 10 _1 M NaCl solution). 



relationships, the low-frequency and high-frequency sur- 
face conductivity terms are expressed as 



-r=^ j4 % ml [%)(i-/)+^]e s „ 

and the parameter A((f>, m) is given by 
1 



(65) 
(66) 

(67) 



Note that the ratio of the surface conductivity to the pore 
water conductivity is often called the Dukhin number (Du) 



in the literature. The expression for A(4>,m) contrasts with 
the high salinity conductivity model developed recently by 
Revil [2012], which yields, 



A{<j>,m) =m(F- 1) 



3 VI -. 



(68) 



For spherical grains, Ftp 3/2 (more precisely 7r/2 for the 
ions to move around a perfect spherical grain) and therefore 
the presently derived expression of A{<j>, m) is 



A(<t>,m) 



\{F-r 



i -, 



(69) 



The difference between equations (67) and (68) is therefore 
really in the factor m in the expression derived from the 
effective medium theory. Note for spherical grain 
(F(f)tt 3/2) and at high porosity, the formation factor is 
given by 



3 (\ 



(70) 



Equations (69) and (70) yields 



A((j>,m)tt\. (71) 
With this approximation and using p = n — 1, the high- 
frequency and low-frequency conductivities are given exactly 
by the Waxman and Smits [1968] equation (see also Waxman 
and Thomas [1974] for the quadrature conductivity) : 



1 n 



°w + /?(+)(! ~f M 



0 s , 



(72) 



*oc = ^s n w {a w + [/3 (+) (l -f u ) + Pf + / M ] J}. (73) 

The present analysis therefore provides a theoretical foun- 
dation to the Waxman and Smits [1968] model in which 



Table 1. Properties of the Berea Sandstone (Clay Type: Kaolinite 


and Illite) 




Property 


Value 


Archie exponent m (-) 


1.78, a 1.69, b 1.85 ± 0.05 e 


Saturation exponent n (-) 


1.83, a 1.98 ± 0.23, e 2.11 -2.17 s 


Specific surface area S sp (m 2 g~') 


0.74 c - 0.93 d 


Porosity (f> (-) 


0.20 ± 0.03 c 


Formation factor F (-) 


18to26 b 


Grain density p s (kg m~ 3 ) 


2660' 


CEC (C kg" 3 ) 


256 C kg" 111 


Excess charge density Q v (C m~ 3 ) 


2.7 x 10 6j - 5.8 x 10 6i 



a Jun-Zhiand Lile [1990]. 

b From the raw data of Lesmes and Frye [2001], corrected for surface 
conductivity. 

"Lesmes and Frye [2001] (one sample). 

d Zhan etal. [2010]. 

"Attia [2005] (eight samples). 

f Aggouretal. [1994]. 

g Z/i» and Toksdz [2012] 

h Using CEC = Q s S sp with Q s = 0.32 C m~ 2 and S sp = 0.8 m 2 g _1 . 
'Vinegar and Waxman [1984, Tables 7 and 10] (porosity 0.19, four 
samples). 
J Usinge F =p s [(l -</»)/0]CEC. 
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Table 2. Properties of the Berea Sandstone 





Permeability 


Pore Throat 


Porosity = <f> (-) 


ko (mD) 


Size 3 A (in urn) (1) 


0.206 b 


498 b 


7.9 


0.205 b 


496 b 


7.9 


0.207 b 


477 b 


7.7 


0.180° 


228 c 


6.0 


0.230 d 


450 d 


6.8 



a From A = 



"k 0 with m= 1.74. 



From Aggour et al. [1994]. 
"From Lesmes and Frye [2001]. 
d From Zhu and Toksdz [2012]. 



their mobility B is related to the mobility j3t + ) by B = 

/%)(! -/m). 

[20] In the high-porosity limit discussed previously, the 
normalized chargeability is given as 



M„ 



(74) 



which is exactly equation (43) of Revil [2012] for the quad- 
rature conductivity obtained through the differential effec- 
tive medium approach. Up to this point, I have not 
discussed the frequency dependence of the electrical con- 
ductivity. This will be done in the next section. Also, rather 
than writing the surface conductivity and the normalized 
chargeability as a function of the CEC, we can use the rela- 
tionship developed by Revil [2012] to express the parame- 
ters as a function of the specific surface area. The 
relationship between the CEC (in C kg -1 ) and the specific 
surface area S Sp (in m 2 kg -1 ) is CEC = QsSs P , where Qs 
denotes the mean charge density for clay minerals (typi- 
cally, two elementary charges per nm 2 , i.e., Q s = 0.32 C 
m , see Revil [2012]). Note that Q s and Q d are related to 
each other by (1 -f M )Qs = Qd- 

2.3. Frequency Dependence of the Effective 
Conductivity and Permittivity 

[21] Figure 3 shows the effective quadrature conductivity 
spectrum for a Berea sandstone (data from Lesmes and 
Frye [2001]), which is a clayey sandstone with a very small 
amount of clay minerals. The properties of the Berea sand- 
stone are summarized in Tables 1 and 2. I call Type A, 
clayey materials showing the type of quadrature conductivity 
spectra shown in Figure 3a. This type of quadrature conduc- 
tivity distribution exhibits three distinct regions : (i) Domain 
I is defined by frequencies such as/ <f P (~ 0.03 Hz for the 
Berea sandstone data shown in Figure 3a). In this region, the 
effective quadrature conductivity decreases with the fre- 
quency (with <j"(uj = 0) = 0). Domain II is defined byf p < 
f <fd (~ 0-3 MHz in Figure 3a). In this domain, the quadra- 
ture conductivity reaches a plateau. Domain III is defined by 
/ >fd- In this domain, the effective quadrature conductivity 
increases with the frequency and is controlled by the permit- 
tivity of the material. I first provide some expressions for the 
two critical frequencies and then I will provide expressions 
for each domain. 

[22] In Appendix B, I show that a characteristic relaxa- 
tion time of a clayey material can be determined by the fol- 
lowing equations 



A 2 

2/>s 

ZLr (+) 



4kF 



v 2 



(75) 



(76) 



where D? + , 
.2 



denotes the diffusion coefficient of the counterions 
(in m z s _1 ), k denotes the permeability (in m 2 ), and A is a 
characteristic pore throat size that can be derived from the per- 
meability and the formation factor (see for instance, Table 2 
for the Berea sandstone and a definition of this parameter in 
Appendix B). The diffusion coefficient is related to the mobil- 
ity of the counterions in the Stern layer, , by the Nernst- 



Einstein relationship Df 



(+) 



k b Tpf +) /\q (+) \ 



where k h is the 



a. Bentonite sand mixture (Slater and Lesmes [2002]) 
3.0 




o 

z 
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- Quadrature conductivity, a" (|j,S/cm) 

b. Saprolite (Revil et al. [2012b]) 
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Quadrature conductivity, a" (uS/cm) 



Figure 4. Relationship between the normalized charge- 
ability and the quadrature conductivity, (a) Eight bentonite 
sand mixtures (R = 0.97). Data from Slater and Lesmes 
[2002]. (b) One saprolite core sample (SI 6, see Revil et al. 
[2012b]). The normalized chargeability is determined from 
the high and low asymptotic values of the (in-phase) electri- 
cal conductivity. The discrepancy at low quadrature conduc- 
tivity may be due to the higher uncertainty in extrapolating 
the high-frequency and low-frequency limits for the electri- 
cal conductivity. 
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Table 3. Properties of the Sandstones From the Vinegar and Waxman [1984] Database Analyzed With the Linear Conductivity Model 
(Shaly Sands, NaCl) 



Sample 


Mineralogy 3 


F (-)" 


ir/OO-'Srrr 1 ) 


<K-f 


CEC d (Ckg- 1 ) 


log fi p e (Cm- 3 ) 


f(~) 


3477 


K, C 


15.1 


30.6 


0.201 


237.74 


6.40 


0.85 


3336A 


K 


23.0 


23.6 


0.205 


393.65 


6.61 


0.93 


3478 


K, C 


18.9 


21.9 


0.181 


417.70 


6.70 


0.94 


101 


K, I 


16.4 


43.0 


0.210 


531.40 


6.72 


0.91 


102 


K, I 


16.2 


13.1 


0.205 


599.84 


6.79 


0.97 


CZ10 


C 


17.8 


72.6 


0.233 


772.91 


6.83 


0.89 


3833A 


C 


13.9 


60.7 


0.224 


1154.1 


7.03 


0.94 


3126B 


I, S 


11.8 


96.2 


0.249 


1446.1 


7.06 


0.92 


3 847 A 


I, c 


47.7 


29.6 


0.126 


754.56 


7.14 


0.96 


3283A 


K, C, I 


19.4 


44.0 


0.186 


1245.8 


7.16 


0.96 


3885B 


I,C 


28.9 


50.4 


0.204 


1676.7 


7.24 


0.97 


3972E 


S, M,I 


21.7 


44.8 


0.187 


1546.7 


7.25 


0.97 


3258A 


I 


49.9 


86.0 


0.128 


2022.1 


7.56 


0.95 


3891A 


C,I 


39.3 


222.4 


0.192 


3325.2 


7.57 


0.93 


3308A 


S, C 


12.9 


381. 


0.270 


5498.4 


7.59 


0.91 


3323F 


S, C 


14.4 


1065 


0.288 


10145 


7.82 


0.89 


3324A 


S, C 


29.9 


932 


0.220 


6622.6 


7.79 


0.84 


3323E 


s, c 


15.1 


1124 


0.281 


9801.6 


7.82 


0.87 


3324B 


s, c 


29.6 


907 


0.204 


7843.3 


7.91 


0.87 


3306F 


s, c 


47.8 


853 


0.171 


7122.6 


7.96 


0.87 



a From Vinegar and Waxman [1984, Table 1], K: kaolinite, C: chlorite, I: illite, and S: smectite. 
b From the fit of the in-phase conductivity data using the linear conductivity model. 
°From the fit of the in-phase conductivity data using the linear conductivity model. 
d From Vinegar and Waxman [1984]. 

e From Q v and porosity assuming a mass density for the grains of 2650 kg itT 3 . 
f Using/= 1 - Q v oo/Q v . 



Boltzmarm constant (1.3807 x 10 23 J K ! ), Tis the absolute 
temperature (in K), and \q(+)\ is the absolute value of the 
charge of the counterions in the Stern layer. Using the data- 
base of Vinegar and Waxman [1984], Revil [2012] found 
/?f +) (Na + , 25° C) = 1.5 x 10" 10 mW -1 . This yields the 
diffusion coefficient of the counterions Z)^(Na + , 25° C) = 

3.8 x 10~ 12 m 2 s -1 for clays. For fully water-saturated materi- 
als, the lower critical frequency at which the effective quadra- 
ture conductivity becomes constant is given by 



A 2 
1U (+) 

Taking A = 6 um and Z)f + >(Na + , 25°C) = 3.8 x 10 



(77) 



- 12 m 2 



the critical frequency associated with the previous 



relaxation time is 



2-KTp TvA 2 



fp 



SnFk 0 : 



(78) 



(79) 



surface of clay minerals is in average Qs = 0.32 C m~ 
(see Revil [2012]). So the average distance between the 
charged sites is 1.8 nm. Using this distance, equation (78) 
yields a critical frequency of 0.3 MHz. 

[24] The second characteristic frequency comes from the 
transition between the low-frequency polarization of the 
electrical double layer and the high-frequency dielectric 
effect in the overall polarization of the material. It is there- 
fore defined by 



Td : 



s'/a", 



(80) 



where a" and e' denotes the quadrature conductivity of the 
plateau and the high-frequency dielectric constant. From 
equation (53), the high-frequency dielectric constant is 
given (for saturated conditions) by 



1 



[e w + (F-l)e s ]. 



(81) 



Taking s w = 80 eo, £ s = 5.9 eo, and F = 18 for the Berea 
sandstone, we get e' = 10.0 eo- The second critical fre- 
quency is therefore given by 



and I obtain f p = 0.034 Hertz in excellent agreement with 
the data shown in Figure 3 a. Note the measurement of this 
relaxation frequency can be used to estimate therefore the 
permeability of the material using equation (79) (this idea 
will be explored below in section 6). 

[23] I follow here closely Vinegar and Waxman [1984] 
in estimating the highest frequency at which the Stern layer 
polarization can play a role. The shortest possible spatial 
variation between charged and uncharged zones is the aver- 
age distance between the surface sites on the clay minerals. 
I have mentioned that the surface charge density on the 



fd 



1 



2nr d 



a 



(82) 



Using the value of the quadrature conductivity for the 
Berea sandstone, a" = — 1.5xl0~ 4 S m _1 (see Figure 3), 
and e' = 10.0 eo yields fd = 0.3 MHz (see Figure 3a). 
Above this frequency, the polarization response of the ma- 
terial is dominated by the dielectric polarization of the ma- 
terial (see Figure 3a). Note also if the quadrature 
conductivity is on the order of — 1 x 10~ 5 S m _1 , the criti- 
cal frequency f d can be as low as 10 kHz. 
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O Revil et al. [2012b], saprolite (rich in illite) 
• Vinegar and Waxman [ 1 984], shaly sands 




Figure 5. High-frequency volumetric charge density ver- 
sus the total charge density. The high-frequency volumetric 
charge density is determined from electrical conductivity 
measurements at various salinities (from the surface con- 
ductivity and the formation factor), while the total charge 
density is determined from the porosity and the cation 
exchange capacity. 

[25] I now describe the quadrature conductivity response 
below the critical frequency f p for Type A materials, which 
are exhibiting a clear plateau in their polarization spectra. 
This plateau is likely resulting from a broad distribution of 
polarization length scales and therefore from the self-simi- 
larity of the pore space for a broad range of scales (this 
idea was expressed by Vinegar and Waxman [1984] and is 
consistent with the modeling performed by Cosenza et al. 
[2008] and observations by Hyslip and Vallejo [1997]). 

[26] For clay minerals, the kinetics of sorption/desorp- 
tion of cations in the Stern layer is pretty fast (few seconds 
according to Pohlmeier and Ilic [1998]). Therefore at low 
frequencies (typically < 1 Hertz), the kinetics of sorption/ 
desorption may influence the polarization of the electrical 
double layer. Wong [1979] provided a general model 
describing the polarization of metallic particles (grains) 
exchanging charge carriers with the surrounding pore 
water. In his case, the charge carriers are redox active spe- 
cies (electron donors and electron acceptors) and the parti- 
cle conductive to electrons. His model is however general 
to work also with insulating grains coated by an electrical 
double layer. Wong [1979] showed that his rather complex 
model can be represented with a Warburg impedance 
model. The Warburg model implies that at low frequencies 



(f < f p ), the quadrature conductivity decreases as the 
inverse of the root mean square of the frequency. For type 
A materials exhibiting a plateau for the quadrature conduc- 
tivity (see Figure 3a), the normalized chargeability is also 
equal to the quadrature conductivity (M„ = —a"). Note that 
Lesmes and Frye [2001] and Slater and Lesmes [2002] 
attempted a semitheoretical explanation for the equivalence 
between normalized chargeability and imaginary conduc- 
tivity. This equivalence is very well illustrated by the data 
of Slater and Lesmes [2002] and Revil et al. [2012b] (see 
Figure 4). 

[27] We therefore propose the following model in which 
the frequency dependence is explicit : 

o , = cr 0 + Mn(j) ' ; f<f P , (83) 

a f = a 00 ; f>f p , (84) 

a" = -M„(£j ' ■ f<f p (85) 

o" = -M n - f>f p , (86) 



Sample, formation factor, surface conductivity 

♦ #3477, F= 15.1±0.2, as= (SO.eiS^xlcf 4 S/m 

A #3833A, F= 13.8±0.6, as= (60.7±12.3)xlo" t S/m 
+ #3323F,_F= 14.4±0.5,av=(1065±89)xl0" 4 S/m 
■ #3323E, F= 15.1±0.7, Os= (1 124±109)xl(T 4 S/m 

• #3308A, F= 12.9±1.0, os= (381±58)xl0" t S/m 




_1 2 r, 1 ♦ 



-4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 0 

log 10 (Du =o s I aj(unitless) 

Figure 6. Pore water conductivity dependence of the con- 
ductivity of shaly sands (data from Vinegar and Waxman 
[1984]). I use here a collection of five samples with more or 
less the same formation factor but different CEC and therefore 
surface conductivities. According to the model developed in 
the main text, the conductivity data should follow the same 
trend, which is the case. The parameter Du denotes the Dukhin 
number, i.e., the ratio of surface to pore water conductivity. 
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and where cfq and are given by equations (63) and (64), 
respectively, and the normalized chargeability M n by equa- 
tion (74). This (empirical) model shows that the quadrature 
conductivity is going to zero at zero frequency [Fuller and 
Ward, 1970; Vinegar and Waxman, 1984]. Therefore, from 
equations (10) to (15), the effective conductivity and the 
effective dielectric constant are expressed as 

ffeff =a 0 +M n (pj ' ; f<f p , (87) 



ffeff = °"oc; / >fp, (88) 



£eff=e' + ^; f>f P , (90) 

where the high-frequency dielectric constant e' is given by 
equation (53). 

[28] Type B spectra are those showing a clear peak in the 
low-frequency part of the spectrum (see Figure 3b for which, 
however, the data are plotted in a linear scale in quadrature 
conductivity). In this case, we expect a Warburg-type model 
to represent the low-frequency complex conductivity [see 
Wong, 1979]. Because the Warburg model is a special case 
of the Cole-Cole model (Cole-Cole exponent c = 1/2), I pro- 
pose the following Cole-Cole model to represent the Type B 
spectra, 



Eeff = £ + 



Mn_(f 



1/2 



/ < fp, 



(89) 



a* (co) = a a 



M 



1 + (i0JT p ) C 



(91) 



• Berea sandstone. Lesmes and Frye [2001] 1 S/mNaCl 
OBerea sandstone. Garrouch and Sharma [1994] 0.4 S/mNaCl 
+ Sandstone. Poley et al. [1978] 3.3 wt. NaCl 



oc-polarization 



cc-polarization + ^-polarization 




eft 



^(K' = £'/£ 0 ) ) 



_l I I I I I I I I L_ 



K 



-I I I I I L_ 



J l_ 



_i i i i_ 



10 



Log (frequency, Hertz) 

Figure 7. Effective relative permittivity versus frequency for the Berea sandstone (data from Lesmes 
and Frye [2001], pH 8, 10" 1 M NaCl solution, Garrouch and Sharma [1994], 2000 ppm NaCl, and Poley 
et al. [1978], sandstone porosity 0.145, salinity NaCl 3.3% wt.). These data are the same than shown in 
Figure 3a but plotted in terms of the effective permittivity rather than in term of effective quadrature 
conductivity. In the geophysical literature, it is very frequent that the a polarization in Domain II is 
(wrongly) assumed to be a Maxwell- Wagner polarization. 
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M n 
1 + (iu)T„ 



(92) 



where r p = l/(2irf p ) denotes the main relaxation time (in s) 
f p denotes the associated critical frequency, and c denotes 
the Cole-Cole exponent (the Warburg model corresponds 
to the case where c = 0.5, and I expect more generally 
that 0 < c < 0.5), and M denotes the chargeability 
(M = ((Toe — 0o)/coc)- The real and imaginary parts of the 
complex conductivity of the Cole-Cole model are given by 
[Cole and Cole, 1941], 



aoo --M n \\ smh[cln(an>)] 

°° 2 "1 cosh [c In (ojTp) ] + cos [f(l 



a" 



M„cos f (1 



2 1 cosh[cln(wTp)] + sin[f (1 - c)] 



(93) 



(94) 



describes well the conductivity data of shaly sands. I use 
five samples from the Vinegar and Waxman [1984] data- 
base with roughly the same formation factor (in the range 
12.9 to 15.1). According to the model, the conductivity 
data normalized by the conductivity of the pore water 
should follow a unique trend when plotted as a function of 
the Dukhin number defined as the ratio of surface to pore 
water conductivity. This is the case in Figure 6, indicating 
that the linear conductivity model described in the present 
paper works well in predicting the pore water conductivity 
dependence of the electrical conductivity of clayey mate- 
rials. This is probably because of the continuity of the dif- 
fuse layer from grain-to-grain as mentioned in section 2. 

4. Low-Frequency Effective Permittivity 

[32] In this section, I test five predictions of the model 
either for the effective quadrature conductivity or the 



Therefore, from equations (10) to (15), the effective con- 
ductivity is given by <7 e ff = d and the effective dielectric 
constant are expressed as 



£eff 



e' + - 



M„ 



cos[f (1 - c)] 



2u> 1 cosh [c In (ojt p )] + sin[f (1 — c) 



(95) 



[29] As mentioned earlier, I expect to have the Cole- 
Cole exponent c in the range 0 (fiat response) to 0.5 (pure 
Warburg response). In the following, I will focus more on 
Type A spectra [see Vinegar and Waxman, 1984; Slater 
and Lesmes, 2002] rather than on Type B [Scott and 
Barker, 2003 ; Binley et al, 2005]. 

3. Effective Conductivity: Theory Versus Data 

3.1. Determination of the Fraction of Counterions in 
the Stern Layer 

[30] The first test concerns the inference of the fraction of 
counterions in the Stern layer, f M . I use equations (63) and 
(65) with equation (69) for the parameter A(<j>, m) to reinter- 
pret the electrical (in phase) conductivity data from Vinegar 
and Waxman [1984] on shaly sands. The results are reported 
in Table 3 and the values of f M determined from the poros- 
ity, the formation factor, and the surface conductivity are 
shown as a function of the excess of charge per unit pore 
volume in Figure 5. The high-frequency excess charge den- 
sity Q v is estimated from the surface conductivity and the 
formation factor by 



Qv 



Q~sF 

/3 {+) A((f),m) 



(96) 



We see that a good average value for f M is close to 0.9 [see 
also Revil, 2012]. This means that approximately 90% of 
the counterions are, in average, located in the Stern layer in 
agreement with triple layer models (see discussion in Revil 
[2012]). 

3.2. Salinity Dependence of the In-Phase Electrical 
Conductivity 

[31] In Figure 6, I investigate if the linear conductivity 
model described in the present modeling approach 



• Binley et al. (2005) (natural water Na + , Mg 2 + Ca 2+ , 0.1 S/m) 
Eggborough and Hatfield sandstones 

• Slater and Glaser (2003) (NaCl, 0.1 S/m) Sandy sediments 
X Lesmes and Frye (2001) (NaCl, 0.1 S/m) Berea sandstone 
+ Revil et al. (2012) (NaCl, 0.1 S/m) Saprolites 

□ Weller et al. (2011) (NaCl, 0.1 S/m) Sandstones 

9 Revil et al. (unpublished, NaCl, 0. 1 S/m) Clean sandstones 
0 Revil et al. (unpublished, NaCl, 0. 1 S/m) Clayey sandstones 
O Data from Bonier (1992) (NaCl, 0.1 S/m) 

□ Data from Koch et al. (201 1) (NaCl, 0.04-0.06 S/m) 




Specific surface area, S (in m 2 /kg) 

Figure 8. Influence of the specific surface area S$ p upon 
the quadrature conductivity. Trend determined for clayey 
sands from the model developed by Revil [2012] at 0.1 
S m" 1 NaCl. Data from Slater and Glaser [2003], Weller 
et al. [201 1], Borner [1992], Koch et al. [201 1], Revil et al. 
[2012b], and Lesmes and Frye [2001]. Measurements 
reported at 10 Hertz for the clayey materials, at 1.4 Hz for 
the data from Binley et al. [2005]. 



317 



REVIL: ELECTRICAL PROPERTIES OF POROUS MATERIALS 



low-frequency effective dielectric constant (these two pa- 
rameters are related to each other by equation (5) in section 
2.1). Figure 7 shows a typical effective relative permittivity 
spectrum. In this section, we are interested to test the model 
for frequency much below the critical frequency f d . 

4.1. Quadrature Conductivity: Influence of the 
Specific Surface Area 

[33] From equations (74) and (86) and for Type A spec- 
tra, the quadrature conductivity of the plateau at saturation 
of the water phase is given as a function of the CEC or the 
specific surface area by 



o" = -(^ Ps 0{ +) fM^jCEC, (97) 

°" = -(^PsPI + JmQs^S Sp , (98) 

where we have used the relationship between the CEC and 
the specific surface area CEC = QsSs p [see Revil, 2012]. 
For clayey sands, taking /3^ +) (Na + ) = 1.5 x 10" 10 m 2 s" 1 
V" 1 at 25°C,/ = 0.90, Q s = 0.32 C m -2 , and p g = 2650 
kg m~ 3 , 1 obtain a" « -aS Sp with a = 7.6 x 10~ S kg m~ 3 
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Figure 9. Influence of the CEC and specific surface area S$ p upon the effective relative permittivity, (a) 
At low frequencies, the relative effective permittivity is linearly correlated with the CEC in agreement 
with the model developed in the main text, (b) At low frequencies, the relative effective permittivity is also 
linearly correlated with the specific surface area S Sp (in m g -1 ) (data from Knight and Nur [1987], various 
sandstones including three Berea sandstones), (c). Measurements at 100 kHz. The line corresponds to the 
trend predicted by the model. Data from Comparon [2005] (100 kHz, clayrocks) and Knight and Nur 
[1987] (105 kHz, sandstones, corrected for the saturation by dividing with the saturation). 
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[Revil, 2012]. For clean sands and sandstones, if we use 
/3 (+) (Na + ) = 5.2 x 10" 8 m 2 s" 1 V"\/= 0.50, Q s = 0.64 
C m" 2 , p g = 2650 kg m" 3 , 1 obtain a = 2.9x 10" 5 S kg m" 3 
[Revil et ah, 2012a]. These two linear trends are shown in 
Figure 8 together with data from the literature. For the data 
of Binley et al. [2005], I converted the 5 por data into specific 
surface area data using Ss p = (S ?w <p/ p g )/(\ — <f>) with a 
grain mass density of p g = 2650 kg m~ 3 and a mean porosity 
of 0.28. The model agrees with the data especially for clayey 
materials. For nonclayey materials, additional data points are 
needed. This demonstrates also the strong difference of the 
mobility of the counterions in the Stern layer for silica sands 
and for clays. This point may be, in turn, related to the density 
of counterions in the Stern layer of silica and aluminosilicates. 
In the case of silicate, only a small fraction of the active sur- 
face sites are locally associated with the sorption of the coun- 
terions in the Ster layer. For clay minerals, I have mentioned 
above that in average 90% of the charged surface sites are 
associated with counterions present in the Stern layer. 

4.2. Effective Permittivity: Influence of the CEC and 
Surface Area 

[34] At low frequency, we can neglect the influence of 
the true dielectric constant in the effective permittivity. We 
first look for an expression of the effective relative permit- 
tivity as a function of the CEC and specific surface area at 
full saturation. We obtain 



Krff {$w = 1) ~ «' +- 



1 
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(99) 



Krff (*w = 1) 



m s 



CEC. 



(100) 



According to my model, the low-frequency effective rela- 
tive permittivity is therefore proportional to the CEC. This 
is in agreement with the experimental results reported in 
Figures 9a and 9c. The slope of the trend shown in Figure 
9a is 1.19±0.13 kg C _1 . Taking m = 2,f= 10 kHz, p s = 



2650 kg m" 3 , $ 



+). 



1.5 x 10" 10 m 2 V" 1 s _1 ,/ M = 0.9, 



the predicted slope is 0.9 kg s _1 C _1 , therefore, in excellent 
agreement with the observations. Equation (100) can be 
also expressed in terms of the specific surface area, 
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which can be written simply as, 
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where b is a constant independent on the specific surface 
area. In Figure 9c, we report measurements of the relatively 
effective permittivity as a function of the specific surface 
area. A fit of these data with equation (102) yields n' = 
25 ± 3 and b = (18 ± 2) x 10 3 kg s" 1 m -2 . 

4.3. Effective Permittivity: Influence of the Saturation 

[35] The third point to check is the dependence of the low- 
frequency effective relative permittivity with the saturation. 



Low frequency behavior 
a. Quadrature conductivity a. Effective permittivity 




Water saturation s w (-) ^ saturation ^ ( . } 

Figure 10. Saturation dependence of the quadrature conductivity and effective permittivity for the 
Berea sandstone, (a) Quadrature conductivity (data from Vinegar and Waxman [1984], 30 Hz, soltrol, 
25°C). Soltrol is an isoparaffin solvent, (b) Relative effective permittivity of the Berea sandstone at two 
distinct frequencies (10 kHz and 100 kHz) showing that the permittivity follows a power-law function of 
the saturation (data from Garrouch [2001], the Berea sandstone is water-wet and saturated by «-decane 
and 10,000 ppm NaCl brine ^2 S m~ , porosity 0.20). Note the consistency of the value of the p expo- 
nent. The plain lines correspond to the prediction of the model. 
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High frequency behavior 
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Figure 11. Saturation dependence of the effective permit- 
tivity. Relative effective permittivity at very high frequen- 
cies (300 and 500 MHz) for the Sherwood sandstone 
(porosity 0.36, data from West et al. [2003]. The plain line 
corresponds to the prediction of the model. 



We consider frequencies low enough so the true dielectric 
constant can be neglected. The saturation dependence of the 
effective relative permittivity is therefore given by 



Keff (s w ) 



o"{s w 



2tv£ 0 



(103) 



where /is the frequency (in Hertz) and a"(s w = 1) repre- 
sents the value of the quadrature conductivity at saturation. 
We therefore see that the permittivity follows a power-law 
function with the saturation. This behavior is in very good 
agreement with the data of Garrouch [2001] (see Figure 
10b). The saturation exponent p is equal to 1.4, fairly con- 
sistent with p = 1.3 resulting from the low- frequency quad- 
rature conductivity measurements of Vinegar and Waxman 
[1984, Table 8] (see Figure 10a). Vinegar and Waxman 
[1984] reported that the saturation exponent n is in the 
range 2.1-2.2 for the Berea sandstone (see Table 1). This 
would yield p = n — 1 ss 1.1 — 1.2 only slightly below the 
previous estimates (1.3—1.4). 

5. High-Frequency Effective Permittivity 

[36] For frequencies much above the critical frequency 
fd given by equation (82) and much below the relaxation 
frequency of water and bound water, the effective relative 
permittivity is given by 



k' = - [ s y w + (f - ik + (i - s n w )*a] , 



(104) 




Porosity (-) Porosity (-) 

Figure 12. High-frequency relative permittivity of the Fontainebleau sandstone, (a) Permittivity meas- 
urements between 50 MHz and 1 GHz (data from Coutanceau-Monteil [1989]). The lines represent the 
prediction of the model, (b) Formation factor/porosity relationship (Zamora et al., unpublished results). 
The plain lines correspond to the prediction of the model. 
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where k w = 80 and K a = 1 denote the relative permittivity 
of water and air, respectively. I now test the predictions of 
equation (104) with respect to saturation, porosity, and CEC. 

5.1. Influence of Saturation 

[37] I use the data from West et al. [2003] at 300 and 500 
MHz, therefore, below the dielectric relaxation frequency 
of water. A fit of the data (shown in Figure 11) yields m = 
1.7 ± 0.1 and n = 1.8±0.1 (for the same sandstone, Binley 
et al. [2002] obtained from low-frequency electrical con- 
ductivity measurements m = 1.7 and n = 2.0, respectively), 
and a relative permittivity for the solid phase of ks = 5.9 
(in agreement with the range provided by Robinson and 
Friedman [2003] for similar sandstones). 

5.2. Influence of Porosity 

[38] In Figure 12, I test the effect of porosity upon the 
high-frequency relative permittivity. I use the data from 
Coutanceau-Monteil [1989] with the Fontainebleau sand- 
stone (99.98% silica). The measurements were performed 
at high frequencies between 50 MHz and 1 GHz. The lines 
in Figure 12 represent the prediction of the model. The 
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Figure 13. Comparison between the model and the ex- 
perimental data from Coutanceau-Monteil [1989] for lime- 
stones, (a) The samples are saturated by fresh water and 
air. (b) Calibration of Archie's law from the measurement 
of the formation factor and porosity for samples of the 
same formation than used for the relative permittivity (50 
MHz). The plain lines correspond to the prediction of the 
model. 



application of the model to the data yields a cementation 
exponent of 1.5 and a relative permittivity for the solid 
phase of 4.5. This relative permittivity is very close to the 
relative permittivity for Si0 2 (3.9 according to Gray et al. 
[2009] and 4.5 according to Olhoeft [1981]). The cementa- 
tion exponent is also consistent with electrical conductiv- 
ity measurements (Figure 12b). In Figure 13, the 
prediction of equation (104) for the relative permittivity/ 
porosity relationship is tested for carbonate rocks. The 
model is consistent with the data and the value of the ce- 
mentation exponent determined from electrical conductiv- 
ity measurements. 

5.3. Influence of CEC 

[39] According to equation (104), the high-frequency rel- 
ative permittivity should be independent on the cation 
exchange capacity. Figure 14 displays various experimental 
data and seems to show that this is the case. 

6. Relaxation Time, Pore Size, and Permeability 

[40] In Figures 15 and 16, 1 test the relationship between 
the relaxation frequency and the pore throat size given by 
equation (78) for clayey sands and clean sandstones, 
respectively. For the Triassic sandstones reported by Scott 
and Barker [2003], I divided the mercury pore diameter by 
2 to get a mercury pore access and then I divided the result 
by 5.3 to get a value of A. Indeed, the capillary entry pres- 
sure is related to the mercury pore radius r c by p e — 2j/r c , 
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Figure 14. At high frequencies (2 MHz), the effective 
permittivity appears independent (or weakly dependent) on 
the CEC, a result that is also in agreement with the model 
developed in the present paper [Garrouch and Sharma, 
1994]. 
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Figure 15. Relationship between the low-frequency time constant and the mean size of the pore throat 
for clayey materials. The size of the thin section images is 1 mm in x. COx stands for Callovo-Oxfordian. 
The plain line corresponds to the prediction of the model. 



where 7 represents the surface tension and r c represents the 
smallest pore of the set of largest pore percolating through 
the porous material. Katz and Thompson [1987] developed 
a relationship between the permeability and the percolation 
length scale r c using percolation principles: ks = r c 2 / 
(226F). A comparison between this equation and equation 
(Al 1) of Appendix B yields r c w 5.3 A. 

[41] It seems that the proposed relationship works well. 
There is however a need for further investigations to con- 
firm this finding, especially for clayey materials. In Figure 
17, 1 use equation (79) to predict the permeability from the 
low-frequency relation time r p and the formation factor F. 
The results show that the permeability can be pretty well 
predicted inside an order of magnitude from the complex 
conductivity. The data covers 1 1 orders of magnitude. Note 
that the diffusivity of the counterions for clean sands is 
equal to the diffusivity of the counterions in the bulk pore 



water (see discussion in Revil et al. [2012a]) while it is 2 
orders of magnitude smaller for clays [see Revil, 2012]. 

7. Relationship Between m, n, and p 

[42] To reduce the number of input parameters to bridge 
reactive transport modeling and geophysical predictions, I 
look now for some empirical relationships between the 
three exponents used in the model. The cementation expo- 
nent m can be determined from the cation exchange 
capacity according to 

m = m Q + m 1 CEC, (105) 

with m 0 = 1.80 and m, = 4.3 x 10" 5 kg C" 1 and the CEC 
is expressed in C kg -1 [Revil et al., 1998]. We test the two 
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relationships to predict the saturation exponent from the ce- 
mentation exponent. The first one is purely empirical and 
has no theoretical foundations but is broadly used in the lit- 
erature in absence of independent estimates of n : 

n = m. (106) 



This relationship is evaluated in Figure 18. Note that a 
number of published papers, in which the two Archie expo- 
nents m and n are reported, cannot be used to assess equa- 
tion (106). Indeed, in these papers, surface conductivity is 
not taken into account and the reported values of m and n 
may be, therefore, grossly wrong. Figure 18 seems to indi- 
cate that the approximation n = m is good enough to be 
used in reactive transport modeling (at least in a stochastic 
sense with an associated probability density). 

[43] The approximation n — m also provides a simple 
equation to predict the water content from the permittivity. 
Taking indeed m = n and neglecting the contribution from 
the air, the permittivity of the air yields 



(107) 



0 Revil et al. (unpublished, clean sandstones) 
■ Scott and Barker (2003) Triassic sandstones 
□ Koch et al. (20 1 1 ) Clean sands 
■ty Binley et al. (2005) Clean sandstones 
♦ Tong et al. (2006) Sandstones 
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Clean sand 



Pore throat size (in urn) 

Figure 16. Main relaxation time t 0 versus the pore size 
A. The pore size data from this study are median values 
obtained from mercury injection data or permeability data. 
For the data from Koch et al. [201 1], the pore size is deter- 
mined from the median grain size and the formation factor 
using the relationship developed by Revil and Florsch 
[2010]. The mean pore size is either determined from mer- 
cury intrusion porosimetry [Binley et al., 2005] or from 
permeability for the data of Tong et al. [2006]. The plain 
lines correspond to the prediction of the model. 



-fit) 
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(108) 



which is checked in Figure 19 for a bentonite-kaolinite 
mixture. 

[44] We can also look for a way to predict the quadrature 
exponent from the saturation exponent. Following Vinegar 
and Waxman [1984], Revil [2012], and section 2 above, I 
can test the following relationship, 



p = n — 1. 



(109) 



In Figure 20, I plot various data from the literature for 
which p and n were independently determined. It seems 
that the quadrature conductivity exponent p is comprised 
between (« — 1) and (n — 0.5). 

[45] If we account for equations (104) to (106), the reduc- 
tion in the number of the parameters is such that the porosity, 
a pore throat size, and the CEC of the material are the only 
three fundamental textural parameters required to predict the 
complete electrical response of the clayey materials in the fre- 
quency range 1 mHz-1 GHz. Of course, the environmental pa- 
rameters (including the salinity and the pH of the pore water 
and the temperature) also control some of the parameters used 
in my model. However, these dependencies are well 
established. 
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• Revil et al. [2012] Saprolite (sample SI 6) 
x Lesmes and Fry [200i] Berea sandstone 
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Figure 17. Predicted versus measured permeability for 
various clean sands and clayey sands and sandstones using 
the low-frequency relaxation time r p or the Cole-Cole 
relaxation time. 
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X Sands, Martinez et al. [2012] 

+ Shaly sands, Clavier et al. [1984] 

■ Sherwood sandstone, West et al. [2003] 



Berea sandstone, Atthia [2005] 
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Figure 18. Saturation versus cementation exponent (n 
versus m). Data from Clavier et al. [1984], West et al. 
[2003], Atthia [2005], and Martinez et al. [2012]. All these 
data are corrected for the effect of surface conductivity 
when such a correction is required. 
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■ Air, clayrock (Comparon [2005]) 
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Figure 20. Saturation exponent p of the quadrature con- 
ductivity/effective relative permittivity versus the second 
Archie exponent n for the in-phase conductivity. The sand 
data are from Schmutz et al. [2012], Vinegar and Waxman 
[1984], Ulrich and Slater [2004], and Abdel Aal et al. 
[2006]. Other data are from Garrouch [2001] and Com- 
paron [2005]. 
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Figure 19. Comparison between the simplified permittivity and conductivity models (expressed as a 
function of the water content, see equations (107) and (108)) and experimental data from Comparon 
[2005] at 1.3 GHz for mixtures of MX80(bentonite)/kaolinite. Porosity 0.40. 
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8. Concluding Statements 

[46] A simple model has been developed to predict the 
apparent conductivity and permittivity (or alternatively the 
apparent quadrature conductivity) as a function of the fre- 
quency, clay content, and clay mineralogy. This model has 
been tested on a broad number of experimental data and 
seems to explain consistently the available data pretty well. 
It also offers the possibility to predict permeability from the 
critical frequency observed at low frequencies in the quadra- 
ture conductivity. My goal with this model is to start to have 
a unified petrophysical framework that can be used to 
develop time-lapse joint inversion algorithms for DC resis- 
tivity, induced polarization, induction-based EM, seismoe- 
lectricity, and GPR for shallow subsurface applications, 
especially regarding the vadose zone and agriculture. 

Appendix A: Saturation Dependence with the 
DEM Approach 

[47] The differential effective medium (DEM) approach 
treats the fluid and solid phases in an asymmetric way : one 
of the two phases has to be the host for the other (for 
instance, grains immersed in a background electrolyte). 
This approach therefore offers a view that complements the 
volume-averaging approach discussed in the main text. The 
solution of the differential effective medium scheme yields 
the following expression for the electrical conductivity 
[Senetal, 1981]: 



1 - g v/°"w 
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(where we have replaced F of the saturated case by Fs~"). 
Equation (Al) has, at low and high frequencies, the follow- 
ing closed- form solutions [Revil et al., 1998] 
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Dukhin number is defined by [Revil et al., 
,2012] 
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[49] Therefore, using equations (65) and (66) of the main 
text (with A(4>, m) = 1), we have 
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[50] An equation similar to equation (Al) can be found 
for the permittivity, 



(A7) 



where Kf denotes the permittivity of the pore fluid. An 
approximation of this equation is 
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where 6 = ks/k/. Now, for unsaturated materials we can 
used the following relationship, 

«/ =*>w+ (1 -s n w )n a . (A9) 
[51] Therefore, the saturation dependence of 0 is 

us 
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Appendix B: Relaxation Time Versus Saturation 

[52] We define in this section the two textural parameters 
F and A. The following canonical boundary value problem 
for the normalized potential T for a cylindrical REV of po- 
rous material of length L can be written as (see Pride 
[1994]) 



V 2 r = 0 in V„ 



n • Vr = 0 in S 



L atz : 
0 atz = 



0 



on S, 



(Bl) 
(B2) 

(B3) 



where z denotes the distance along the axis of the cylindri- 
cal core. In these equations, n denotes the unit vector 
normal to the pore water/mineral interface S and V p corre- 
sponds to the pore volume. The boundary conditions defin- 
ing the normalized potential T are representative for the 
injection of an electrical current into a rock sample in the 
absence of surface conduction along the pore/water inter- 
face [see Johnson et al., 1986; Avellaneda and Torquato, 
1991; Pride, 1994]. The dynamic pore radius A and the 
formation factor F are defined as [Johnson et al, 1986] 
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where V is the total volume of the considered REV. The 
length scale A therefore corresponds to a weighted version of 
the hydraulic pore radius VpIS weighted by the norm of the 
electrical field (normalized by the electrical field imposed 
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from the boundary of the material) and therefore A can be 
seen as the (dynamic) characteristic pore size of the material. 

[53] Using a volume-averaging approach, Johnson et al. 
[1986] and Pride [1994] developed the following general 
equation for the electrical conductivity of a water-saturated 
rock: 



00 



(B6) 



where A is a characteristic pore length scale defined above 
and 'Ed denotes the specific surface conductivity associated 
with the diffuse part of the electrical double layer (in S). 
Similarly the high-frequency conductivity and the normal- 
ized chargeability should be obtained by 



1/2 



(B7) 
(B8) 



where denotes the specific surface conductivity associ- 
ated with the Stern layer (in S). To be compatible with 
equations (72) to (74), the following scaling laws for the 
dependence of the formation factor and length scale A with 
the relative water saturation should hold, 



F 
A 



Fs: 



As„ 



(B9) 
(BIO) 



[54] The quasi-static permeability k (for oj = 0) is related 
to the formation factor F and the dynamic pore radius A as 
[Johnson et al., 1986], 



A* 
8F' 



(Bll) 



[55] Therefore, the permeability should scale with the 
water saturation according to the following power law, 



A 2 
SF" 



„2+n 



(B12) 



[56] According to this scaling, the permeability can be 
computed as the product of a permeability at saturation k s 
and a relative permeability that depends only on the rela- 
tive water saturation k(s w ) = ksk r (s w ) with k r (s w ) = s 2 ^" . 
This type of power-law behavior is also the one predicted 
by the Brooks and Corey [1964] model. 

[57] The relaxation time for the low-frequency polariza- 
tion is given by [Revil et al, 2012a] 



A 2 
lu X+) 



(B13) 



[58] Consequently, using equation (BIO), this yields the 
following dependence of the relaxation time with saturation, 



A 2 ^ 



(B14) 



[59] This equation provides a way to compute the low- 
frequency relaxation time from the pore size of the throat 
and from the saturation. 
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